The interplay between rf-wave superconductivity and antiferromagnetic fluctuations: a 

quantum Monte Carlo study 



F.F. Assaad 
Institut fiir Theoretische Physik III, 
Umversitdt Stuttgart, Pfaffenwaldring 57, D-70550 Stuttgart, Germany. 



oo 
in 



CD 



I 

o 
o 



> 

o 

m 

o 

00 



o 



X 



We consider the repulsive Hubbard model on a square lat- 
tice with an additional term, W, which depends upon the 
square of a single-particle nearest-neighbor hopping. At half- 
band filling, constant W, we show that enhancing U/t drives 
the system from a d-wave superconductor to an antiferro- 
magnetic Mott insulator. At zero temperature in the super- 
conducting phase, spin-spin correlations follow a powerlaw: 
g-f Q|^-a jjg].g Q — (tt^tt) and a is in the range 1 < a < 2 
and depends upon the coupling constants W and U. This 
results is reached on the basis of large scale quantum Monte- 
Carlo simulations on lattices up to 24 x 24, and is shown to 
be independent on the choice of the boundary conditions. We 
define a pairing (magnetic) scale by the temperature below 
which shoH range d-wave pairing correlations (antiferromag- 
netic fluctuations) start growing. With flnite temperature 
quantum Monte Carlo simulations, we demonstrate that both 
scales are identical over a large energy range. Those results 
show the extreme compatibility and close interplay of antifer- 
romagnetic fluctuations and d-wave superconductivity. 



The understanding of the interplay between d-wave su- 
perconductivity and antiferromagnetism is a central issue 
for the understanding of the phase diagram of High- Tc 
superconductors 0| . The aim of this work is to further 
study a model which shows a quantum transition between 
a d-wave superconductor and an antiferromagnetic Mott 
insulator. It thus enables us to address the above ques- 
tion. The model we consider, has been introduced in Ref. 
a. It is defined by: 



with the hopping kinetic energy 
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Here, W > 0, <5 = ±c 



, ±a„, and nr = ct cr where 



cl (cj^) creates (annihilates) an electron on site i with 
z-component of spin a. We impose twisted boundary 
conditions: 



(3) 



with $0 = he/ e the flux quanta and L the linear length 
of the square lattice. The boundary conditions given by 
Eq. (^ account for a magnetic flux threading a torus on 
which the lattice is wrapped. One major advantage of the 
above model consists in the fact that at half-band flUing, 
the sign problem in the quantum Monte-Carlo (QMC) 
method may be avoided. This statement is valid even for 
the above boundary conditions. Hereafter, we will only 
consider the half-flUed case. 

Previously we have considered the model of Eq. (|^) at 
fixed values of U and as a function of W. At U/t = 4 
we have shown that the quantum transition between 
the Mott insulator and c?-wave superconductor occurs at 
Wc/t 0.3 We have equally considered the model at 
finite doping and given numerical evidence of the occur- 
rence of a doping induced quantum transition between 
the Mott insulator and c?-wave superconductor Here, 
we fix W/t = 0.35 and vary U/t for the half-filled case. 
The advantage of this choice of parameters is that it 
provides us with a large parameter range in the super- 
conducting phase where QMC simulations are extremely 
precise, thus allowing us to reach large lattices. This al- 
lows us to reliably study the nature of the spin degrees 
of freedom in the superconducting state. We have used 
two QMC algorithms, i) The Projector QMC algorithm 
which produces zero temperature results in a canonical 
ensemble and ii) the finite temperature grand canon- 
ical QMC method 1^,0. The application of those algo- 
rithms for the above model has been discussed in refer- 
ence Both algorithms, generate a systematic error 
proportional to (Ar)^ where At denotes the imaginary 
time step. To determine the exponent of the spin-spin 
correlations we have extrapolated At to zero. 

We first concentrate of the charge degrees of freedom 
at zero temperature. To distinguish between a super- 
conductor and insulator, we compute the ground state 
energy as a function of the twist in the boundary con- 
dition: Eq{^). For an insulator, the wave function is 
localized and hence, an exponential decay of AEo{^) = 
i?o($) — Eq{^o/2) as a function of lattice size is expected 
Q. In the Hartree-Fock spin density wave (SDW) ap- 
proximation for the half-filled Hubbard model, one ob- 
tains Ai?o(<&) = Q!($)Lexp (— L/^) where ^ is the lo- 
calization length of the wave function. On the other 
hand, for a superconductor, Ai?o(*&) shows anomalous 
flux quantization: Ai?o($) is a periodic function of <& 
with period $0/2 and a non vanishing energy barrier 
is to be found between the flux minima |-|l|] so that 
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Ai?o(<I'o/4) remains finite as L ^ cxd. Fig. ([l]) plots 
AEo{^/4) versus 1/L at W/t — 0.35 and for various val- 
ues of U/t. At values oiU/t < 4, the data is consistent 
with a, 1/L form and scales to a finite value. At U/t = 5, 
the data may be fitted to the above SDW form. Thus, 
for values of W/t — 0.35, the data is consistent with the 
occurrence of a quantum transition between a supercon- 
ductor and a Mott insulator at 4 < Uc/t < 5. 

In order to determine the symmetry of the supercon- 
ducting state, we have computed equal time pairing cor- 
relations in the extended s-wave and d„2_,.2 channels: 



with 
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Here, f,{5) = 1 and fd{6) - 1(-1) for S - iai i±ay). 
The vertex contribution to the above quantity is given 
by: 



Plsir^ = PdAr^ - E fdAS)fdAS') (6) 

<j,S,S' 

Per definition, g{r) = in the absence of interactions. 
Fig. ^ plots P^ir) along the diagonal of the lattice. We 
consider lattices ranging up to L = 24. For a fixed lat- 
tice size, one notices a plateau structure as a function of 
distance. The extrapolation of this plateau value to the 
thermodynamic limit is hard and the above criterion of 
fiux quantization, proves to be a more efficient method 
to conclude superconductivity |T2[ . In comparison to the 
d-wave signal, the extended s-wave signal at large dis- 
tances (data not shown) may not be distinguished from 
zero within our accuracy. The extended s-wave pair-field 
correlation dominate at short distances. The data con- 
firming this statement may be found in the table. 

The added PF-term, contains no processes which ex- 
plicitly favor d-wave superconductivity. On the contrary, 
one would expect the W^-term to favor extended s-wave 
symmetry, and this shows up on short length scales. The 
fact that d-wave symmetry dominates at long-range is a 
result of the underlying magnetic structure. 

To study the spin degrees of freedom, we have com- 
puted equal time spin-spin correlations 



S{r) = tM^W)) 
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where S{r) denotes the spin operator on site r. Fig. |^, 
plots S{L/2, L/2) versus 1/L where L corresponds to the 
linear size of the lattice. We consider periodic boundary 
conditions, and various values of U/t. W/t is constant 



and set to W/t = 0.35. For values of U/t < 4 the data is 
consistent with a powerlaw decay: 
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At U/t — 5 extrapolation of the data leads to a finite 
staggered moment, and thus the presence of antiferro- 
magnetic long-range order. The statement in Eq. (^) is 
surprising and deserves confirmation. To cross-check the 
validity of the above equation, we demonstrate numeri- 
cally that it is independent on the choice of the bound- 
ary condition. Fig. | plots S{L/2,L/2) at U/t = 1 and 
U/t = 2 for W/t = 0.35. At $ = 0, the solid lines 
corresponds to least square fits to the form cL~°'. We 
thus determine the exponent a. The data at $ = <I'o/2 
is consistent at large distances with the form cL^°' thus 
showing that the exponent a is independent on the choice 
of the boundary. For the $ = $0/2 simulations we were 
able to reach lattices up to 24 x 24. At U/t = 2 and for 
the boundary conditions set by <i> = $o/2, P^(L/2,L/2) 
(see Fig. |) decays slower than S{L/2,L/2) (see Fig. |) 
thus confirming that the pairing correlations are domi- 
nant [|l3|. 

Having put Eq. (^) on a numerically firm basis, we now 
argue why it is a surprising result. In two dimensions, 
correlation functions which decay slower than l/r^ lead 
to divergences in Fourier space. The Fourier transform 
of the spin-spin correlations a,t U/t — 2 is plotted in Fig. 
||. One sees a systematic increase of S{Q = (tTjTt)) as a 
function of system size [Q. In contrast, a mean field d- 
wave BCS calculation, yields spin-spin correlations which 
decays as a powerlaw with ascs = 3-5. This mean- 
field result leads to a finite S{Q) in the thermodynamic 
limit. Another surprise comes from the dependence of the 
exponent a on the coupling constants, U and W. This 
is a feature which occurs in one dimensional quantum 
systems such as the t — J or Hubbard models. 

We now consider the model at finite temperatures with 
the use of the grand-canonical QMC algorithm. In Fig. ^ 
(a) we plot S(Q) for values of U/t ranging from U/t = 
to U/t = 8. We define the magnetic scale, Tj, as the 
temperature scale where S{Q) starts growing. In Figs. 
I (b) and (c) we consider the vertex contribution to the 
d-wave paring correlations at distance R = {L/2, L/2): 
PAL/2, L/2). Here, we consider an L = 10 lattice. At 
this distance the s-wave pairing correlations are negli- 
gible. We define the d-wave pairing scale, T^, by the 
temperature scale at which PAL/2, L/2) starts growing. 
From the data, we conclude that 



(9) 



That is, antiferromagnetic fluctuations as well as d-wave 
pairing fluctuations occur hand in hand at the same en- 
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ergy scale. If U/t > Uc/t ( U/t < Uc/t ) antiferro- 
magnetic correlations (d-wave pair-field correlations) will 
dominate at low temperatures. We note that in the large 
U/t limit, Tj should scale as t'^/U. That up to U/t = 8 
we still do not see this behavior, may be traced back to 
the fact that the W-tevm enhances the band- width [ p5[ . 
Tp = Tj is a natural consequence of the assumption of 
iS'0(5) symmetry which unifies antiferromagnetism with 
d-wave superconductivity The data, however, does 
not demonstrate the presence of this symmetry in this 
model and further work is required. 

In conclusion, we have considered aspects of the t — 
U ~ W model defined in Eq.(|l|). The model, at half- 
band filling, shows rich physics and allows us to study 
the interplay between magnetism and c?-wave supercon- 
ductivity. We have pointed out the surprising nature of 
the spin-spin correlations in the superconducting state: 
S{r) ^ e^'^ '5|r|^", where Q — (7r,7r), a is in the range 
1 < a < 2 and depends on the coupling constants W 
and U . Those conclusions are based on large scale calcu- 
lations for system sizes up to 24 x 24. We have equally 
shown that the energy scales at which d-wave pairing and 
antiferromagnetic fluctuations occur are identical. The 
further understanding of spin and charge dynamics, as 
well as the doping of the model remains for further stud- 
ies. 
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FIG. 1. A£'o($o/4) = _Bo($o/4) - £o($o/2) as a function 
of inverse linear size L. For U/t < 4 the data is a consistent 
with a 1/L form which extrapolates to a non- vanishing value 
in the thermodynamic limit. At U/t = 5 the data may be 
fitted the SDW from L exp(— L/^) thus signaling an insulating 
state. 
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FIG. 2. Vertex contribution to the d-wave pair field corre- 
lations. Here, the temperature is set to T = 0. 
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W/i = 0.35, {n) = 1 




L = 10, ir/i = 0..33, {n) = 1 
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FIG. 3. S{L/2, L/2) versus 1/L. The simulations were car- 
ried out at ^/^o — 0. Here the At systematic error is ex- 
trapolated to zero and the temperature is set to T = 0. The 
solid lines are least square fits to the form: L~" . The result 
of the fit is summarized in Eq. (H). 
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FIG. 4. S{L/2, L/2) versus 1/L for different boundary con- 
ditions. The solid line for the $ = data is obtained from 
least square fit to the form L~" . The same value of the ex- 
ponent a was used to fit the data at $ = $o/2. Here the At 
systematic error is extrapolated to zero, and T — 0. 
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FIG. 5. Fourier transform of the equal time spin-spin cor- 
relations. 
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FIG. 6. (a) S{Q — (vr,7r)) as a function of temperature, for 
different values of U/t. (b)-(c): P^{L/2, L/2) versus temper- 
ature. The calculations presented in the figure were carried 
out at $ = 0. 



r 






(0,0) 
(0,1) 
(0,2) 


0.2950 ± 0.0018 
0.0932 ± 0.0009 
0.0076 ± 0.0002 


0.1304 ±0.0011 
0.0238 ± 0.0006 
0.0252 ± 0.0003 



TABLE I. Short range vertex contribution of pair-field 
correlations in the extended s- and d-wave channels. Here 
we consider an L = 24 lattice at W/t = 0.35, U/t = 2 and 
(n) = 1. The boundary conditions are set by <1> = $o/2. The 
distance r is in units of the lattice constant. 
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